AXIAL VECTOR CURRENT MATRIX ELEMENTS AND QCD SUM RULES 



m 
o 
o 

(N 
o 

Q 



> 

O 

cn 

m 
o 

i!^ ■ 

Oh, 
I 

Oh! 



X 
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The matrix element of the isoscalar axial vector current, 117^7511 + 75 d, between nucleon states 
is computed using the external field QCD sum rule method. The external field induced correlator, 
(0|g7p75g|0), is calculated from the spectrum of the isoscalar axial vector meson states. Since it is 
difficult to ascertain, from QCD sum rule for hyperons, the accuracy of validity of flavour SU(3) 
symmetry in hyperon decays when strange quark mass is taken into account, we rely on the empirical 
validity of Cabbibo theory to dertermine the matrix element -117^7511 + d^^^sd — 2s7^75S between 
nucleon states. Combining with our calculation of ti7p75M + d^^ysd and the well known nucleon 
/3-decay constant allows us to determine {p,s\ ^uj^j^u + ^djfj.jsd + |s7ai75S \p, s) occuring in the 
Bjorken sum rule. The result is in reasonable agreement with experiment. We also discuss the role 
of the anomaly in maintaining flavour symmetry and validity of OZI rule. 

PACS numbers: 



I. INTRODUCTION 

The determination of flavour singlet and non-singlet 
axial vector matrix elements between nucleon states is 
of considerable interest theoretically and experimentally. 
While non-singlet currents appear in nucleon and hy- 
peron beta decays, the linear combination 



(p, s\ -u-i^,-i^u + -^7^75^ + 7757^755 b, s) ^ -s^,Gbj 



(1) 

appears in the integral of the first moment of the po- 
larised structure function 51(2;) in the sumrule of Bjorken 
Pi . Here \p, s) denotes proton state of momentum p and 
polarisation s^. We have introduced the constant Gsj 
to denote the linear combination of axial vector cur- 
rent with coefficients equal to square of quark charges. 
Nearly three decades ago Ellis and Jaffe using a sim- 
ple minded application of the OZI rule, set the strange 
quark current matrix element between proton states to 
be zero, which immediately enabled them to write 



1 5 
•^Bj = -^Ga + —Gs 



(2) 



where, Ga is the isovector axial vector coupling occuring 
in nucleon beta decay and Gg is the octet current cou- 
pling. Eqn.(|2Jl is in disagreement with experiment if 
SU(3) flavor symmetry is a good symmetry Already 
in 1979 Gross, Trieman and Wilczek had pointed out 
that because the light quark masses are unequal, i.e. 



nid - rriu 



md 



Oil) 



one should expect large violations of Isospin symmetry in 
the Bjorken sumrule if anomaly is neglected. The matrix 
elements of the anomaly between vacuum and Goldstone 
states Itt) and jry) are not zero and play a crucial role 
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in maintaining flavour symmetry. As we shall discuss 
further in Section^ this at once implies that Ellis- Jaffe 
assumptions of simulataneous validity of SU(3) flavour 
symmetry and OZI rule are mutually incompatible. The 
octet current has no anomaly while the singlet does. 

The computation of matrix elements of the axial cur- 
rent in Eqn.([3l is clearly a problem in QCD. This can be 
addressed using the external field method introduced in 
the context of QCD sumrules by loffe and Smilga and 
Balitsky and Yung 7] in 1983 to calculate the magnetic 
moments and which has since then been used for numer- 
ous other matrix elements as well. In fact computation of 
the axial vector current matrix elements has been done 
using these methods in ref.|3|-p^. 

The reasons for reconsidering the earlier QCD sum rule 
determination of Gbj are as follows. 

1. Besides the usual Lorentz invariant chiral and gluon 
condensates present in QCD vacuum, additional 
non- invariant correlators {0\qj^i"f5<l\0) induced by 
the external field, enter the sumrules. If in Eqn.(l), 
on the left hand side, we had only flavour non- 
singlet currents, the corresponding external field 
induced correlator of dimension three, 

(0|m7^75U - J7^75d|0)|Ext. Field and (3) 

(0|m7^75U -I- ^7^75^ - 2s7^75s|0) Isxt. Field (4) 

for example, can be found using the ward identity 
and PCAC. 

The determination of either the isoscalar current or 
SU(3) singlet current correlators 

(0|U7^75M -I- d7^75rf|0) Isxt. Field (5) 

(0|u7^75M -I- djf^j^d + s7^75s|0)|Ext. Field (6) 

with the help of Ward identities are no longer sim- 
ple since they involve the gluon anomaly. In ear- 
lier works some authors ^3, 0, used simply 
the values of the non-singlet current induced cor- 
relators, as an expedient measure, while loffe and 
Khodjamirian [l^ found inconsistencies with SU(3) 



2 



flavor symmetry in their calculations of the correla- 
tor in Eqn.(6). Indeed later in ref.'l5"|, the correla- 
tor Eqn.(6) was treated as an unknown free param- 
eter and its value fixed by using the experimental 
value of Gbj ■ In the present work we follow a dif- 
ferent procedure. The external field induced corre- 
lator corresponding to the isoscalar current Eqn.(5) 
is determined directly from the axial vector meson 
spectrum and we verify that similar determination 
of the nonsinglet isovector current induced corre- 
lator Eqn.(3) indeed yields a value consistent with 
Ward identity and PCAC. 

2. We also clarify the differences in the calculations of 
the Wilson coefficients between ref. 9] and 11] on 
one hand and ref. [T3| on the other. 

3. In the analysis of sumrules, we use a value of the 
QCD scale parameter A consistent with present 
data, while earlier works used a significantly lower 
value. 

This paper is organised as follows. In the next section 
we briefly recapitulate the arguments of Gross, Treiman 
and Wilczek, regarding the role of the anomaly in main- 
taining the flavour symmetry. Although the anomaly is 
superficially a flavor singlet its matrix elements between 
the vacuum and the Goldstone state |7r°) and jyy), are not 
zero. Wc also make a brief digression on the OZI rule. 
In section IIIII we give a summary of the external field 



method, and an appendix explains the differences be- 
tween ref. and [l3| in the computation of the Wilson 
coefficients. In Section Hvl we outline the determination 
of the external field induced vacuum correlators which 
is used in Section to determine the isoscalar matrix 
element and we end with a brief discussion. 



II. FLAVOUR SYMMETRY, OZI, AND THE 
ANOMALY 

We briefly recall the argument of ref . If we ignore 
the anomaly then we have, for the divergence of isoscalar 
current 

-I- z(to„ - md) [U75U - ^75^] 

(7) 

We also have, for the isovector 

Ft,w?^(I)t,q = [1*7^ 75 u - d^filbd] 

= i{mu + md) [u75U - ^75^] 

-I- i(m„ - md) [U75U -I- d^i^d] (8) 

Combining these using PCAC one gets 



(A^|a^(u7^75U-Kd7^75d)|Ar)/=i = 
m„ - md 



m„ -I- md 



(iV|a''(i27^75W-d7^75d)|7V) 



(9) 



where the subscript / = 1 on the nuclcon states denotes 
the difference between proton and neutron matrix ele- 
ments 

{N\0\N)i=, = {p\0\p) - {n\0\n). 

Eqn.® implies a large violation of isospin in Bjorken 
sum-rule since 



md - m„ 
mu + md 



0(1). 



This conclusion is avoided by noting that one has ignored 
the anomaly. In Eqn.(7) one should write 

[u7^75U + d^^_i^zd] = i{mu + m<j) [U75M + ^75^] 
-I- z(to„ - md) [M75M - ^75^] 



(10) 



where GJJ^ — \t^ivafiG°'°'^ ■ Using a Sutherland type ar- 
gument it is derived in ref 



«(0|(m„u75U -I- mrfd75(i)|7r°) 



IGtt 



2 ^ iiV^ iiv\ 



(11) 



and 



^ /^o. r^a 1^0 
lOTT"^ 



(0|2i(m,g75g)|^° 



g^.g;:,!^") = (12) 



where, q = s, c, ... etc. Using again a PCAC argument 
they 5] obtain 

2i(0|m„ M75M -I- J75d|7r°) — — F-,^ m^^V2 

mu + md 

= 4i{0\ms S75s|7r°) 
= 4i(0|TOc C75c|7r°) 
= 4i{0\mb 6756^°) ...(13) 

In other words the matrix elements of the anomaly are 
far from being a flavour singlet. It was pointed out 
by Novikov et.al that the matrix elements of the 
anomaly between vacuum and \r]) is again not zero. 
Writing 

(0|u7p75M + ht^lbd - 2s7^75s|77) = i\/&Fj,k^, (14) 

and taking its diveregence one gets 

(0| - 4zm,s75s|77) = V^F^m^, (15) 

where we have ignored m,i and md in comparison with 
mg. For the singlet current one may write 

(0|"7m75W + d-^^ilbd + s7M75sh) = ^flk^, (16) 

If SU(3) flavour symmetry were exact, as happens when 
all light quark masses are neglected, /i = 0. On the other 
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hand F-^ remains finite in the chiral Uniit of massless 
quarks. We can expect then F^^ » fi. Setting fi = 0, 
it is then easy to obtain from Eqns. (15) and (16) cf.fl^. 



Sets 



(17) 



We learn from Eans. (|13|l and H17|l above that the ma- 
trix elements of the anomaly between the vacuum and 
non-singlet goldston boson |7r°) and j?/) is not zero but 
proportional to quark mass differences. 

We note that lofTe and Shifman [l3 using Eqns. l(T^ 
and H17|l above obtained the result 



r(V>(23) ^ J/7/;(ls)7r°) 
r(V'(2s) ^ J/t^Mr,) 

md-mu\ I mT,\ ( Ptt 



Pr, 



(18) 



Experimentally one has r = (3.07 ±0.70) x 10 ^. We use 
this to find niu/md and obtain 

^ =0.44 ±0.07 

ma 

which is consistent with values obtained using entirely 
different inputs; Gao ,18] rriu/md = 0.44 and Leutwyler 
Ei mu/md = 0.553 ±0.043. 

Encouraged by the agreement between Ean. (|17() and 
experiments let us consider the matrix element 



(0|s7M75sh) = ifsk^. 
Taking the divergence 



(19) 



47r 



Using Eans. p5|l and p7() we have 

-V6 



3 



F^ 



(20) 



(21) 



Now it is known that the Goldberger-Treiman relation 
V2mNGA = g-irNFT, (22) 

is accurate to a few percent and is exact in the chiral limit 
of massless quarks. Corrections to it have the structure 
|2i 



Ciml + 027711171(7111) + .. (23) 



In other words, GT relation is obtained by retaining only 
the Goldstone pole and discarding the continuum contri- 
bution in the dispersion integrands. In the chiral limit, 
equations analogous to Eans. H22|) and (23) also hold good 
for the nucleon matrix element of the octet current with 
Ga replaced hy Gjj + Gd — ^Gs and at replaced by g^jN 



etc. with leading corrections proportional to quark mass 
as in Ean H23(l . If we now naively extend these dispersion 
relation considerations to the nucleon matrix element of 
the strange quark current 57^753, retaining only the 77- 
pole and discard the continuum as well as ?/ pole we 
immediately obtain from Eqn.(14), (19) and H21|l 



or 



Gs — ~iGu 



Gij + Gi 



Gd - 2Gs) 



Gs 







(24) 



(25) 



which is same as the Skyrme model result 21]. 

Our main point in this section is that OZI cannot 
be applied naively. The anomaly is important to avoid 
gross violations of SU(3) flavour symmetry and matrix 
elements of the anomaly are not flavour symmetric. 

It is worth emphasising that OZI rule violates unitar- 
ity, a cardinal property of all S-matrix elements. Correc- 
tions to OZI rule can be estimated u sing unitarity and 
they are very much process dependent |22l ] . Charmonium 
decays illustrate the point. For example. 



B.R.(J/^(ls) 
B.R.(V'(2s) 



Ptt) 
pit) 



= 1.27 X 10"^ 
< 8.3 X 10"^ 



while 



despite the fact Ttot{i^{'^s)) — 277 keV is just a fac- 
tor of three larger than Ttot{J /'^{^s)) — 87 keV. Again 
in the decay of J/ip{ls) into light mesons, SU(3) fla- 
vor symmetry works better in pseudoscalar vector de- 
cays than in vector tensor decays. Also the decay mode 
J/^{ls) (^TTTT is not doubly suppressed as one might 
naively expect. These emphasise the point that correc- 
tions to OZI rule are to be studied individually for each 
matrix element j22| and there is no universal principle 
which tells a priori how good is the OZI rule for any 
specific matrix element. 



III. QCD SUMRULE WITH AN EXTERNAL 
FIELD 

We shall follow closely the notations of loffe 0] . We 
consider the nucleon correlator in an external field 

n(p, A,) =1 I d^x e^P-(OlT{r;(a:), f;(0)}l0) (26) 



where 77(2:) is the nucleon current 

7?(x) = e'^^'^u'^ix) C-i^u\x)n^d-{x) (27) 

with proton quantum numbers; u°',d}' are quark fields 
and a, b, c are color indices. Afj_ refers to constant external 
field. To compute the matrix element of the current 
between a proton state (plj^Jp) one adds a term 



AC = j^A^ 



(28) 
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to the Lagrangian and evaluates 11 (p, A,,) in Ean. l26() 
upto terms linear in A^. In ref the following four 

different currents were considered. 

jl = uii-ilhU - d^^^ad isovector (29) 

j^L = uiiilbU + d-f^-f^d isoscalar (30) 

= u'yti^5U + d'yf,^5d~2sjfj,^5S octet (31) 

jI = -"7^75^ + ^7^75^+57^755 SU(3) singlet(32) 

Bearing in mind that the corrections to OZI rule are a 
priori unknown and can be large, we shall consider only 
the isoscalar current Eqn.(30) here, since the nucleon cur- 
rent of loffe in Eqns. (26-27) above has only up and down 
quark fields, and therefore couples to the 57^755 term 
in the octet current, Eqn.(31), or singlet, Eqn.(32),only 
through giuonic corrections that is by corrections to OZI 
rule. 

In deriving the QCD sumrules one must take into ac- 
count, external field induced correlators . We define 

(0197^75910) U, = FA^ (33) 

{0\gsqYG;^l'l\0)K = HA^ (34) 

It is clear that the constants F and H are in general 
different corresponding to the different AC introduced in 
Eq. (28-32). The constant F can be obtained from 

(0|g7M759|0)U, = ^ J d^x{0\T{AC,q^^Toqm (35) 

and the constant H from 

(0|5.g-yG;^7''9|0)U, 

J d'x{0\TiAC,gsq^G';,^Yqm (36) 

Complete details of the calculation of n(p, Ag) in 
Eqn.lPJ) can be found in ref. 8, 9, KHHIlllil for both 
the isovector and isoscalar currents. For the isoscalar ax- 
ial vector matrix element the sum rule then reads 

bAP 167r2 AP 4 /9 

= ~X%{G + AAP)exp[~m%/M'^] (37) 

Here tun is the nucleon mass and AP is the Borel mass 
variable. In the right hand side, Ajv is defined by 

{OMx)\p) ^ \nv{p) and = ^ 

Gjj + Gd = G is the isoscalar axial vector current matrix 
element. The term AAP arises from the fact that in the 
presence of external field there are non-diagonal transi- 
tion between nucleon and excited states. In the left hand 



side 

ln{AI/A) 2/- \ , / 2^a \ 

Ei^l-il + WyAP) e-W'VM^^ 
= 1 - (1 + W^AP + Wy2AP) e-^'/^^' 

/i is the rcnormalization scale, which we take to be 1 
GeV, and A is the QCD scale which for 3 flavor case is 
247 MeV 113. 

Although the sum rule has been written down earlier 
by several authors the purpose of reconsidering it here 
are the following. 

1. It is clear the constants F and H should be deter- 
mined from other sumrules or different considera- 
tions before we can use it to find G in Eqn.(37). 
As will be explained in detail in the next section, 
for the isovector case, it is relatively easy to deter- 
mine them using PCAC. For the singlet currents 
one must take into account the anomaly. A decade 
ago loffe and Khodjamirian Tsj attempted to com- 
pute F using the Ward identity and sum rules for 
the divergence of the SU(3) singlet current. They 
found inconsistency with SU(3) flavour symmetry 
chiefly because of the large difference between the 
strange quark mass and the up or down quark mass. 
In later works 0, 0| this lead loffe and his col- 
laborators to use the sumrule, Eqn.(37), to find F 
using the experimental value of Gbj in Eqn.(l). 
In contrast here for the isoscalar matrix element 
we shall determine the constant F from the spec- 
trum of isoscalar axial vector mesons, and use that 
value to determine from Eqn.(37) the value of G 
and therefore Gbj- 

2. The external field Lagrangian, modifies the prop- 
agation of the current ri{x), in the QCD vacuum 
in two ways. One in which the external field di- 
rectly couples to the fields in 'ri{x). The other, in 
which the external field modifies the vacuum as 
represented by the induced correlators F and H, 
which appear in the second and fourth terms in 
the sum rule, Eqn.(37). Now consider the differ- 
ence between the three AC in Eqns. (30), (31), and 
(32). As long as gluon loops are not included, the 
additional term 57^755^^ will not couple to the u 
and d fields in i]{x). The strange quark term will 
only affect the constants F and H in Eqn.(35) and 
(36). Putting it differently with this assumption, 
the same sum rule, Eqn.(37), is valid for all the 
three currents, Eqn.(30)-(32) namely, the isoscalar, 
octet and SU(3) singlet with only the external field 
induced correlators F and H being different. Cor- 
respondingly if we were to compute only the ma- 
trix element of (p, s|s7^75s|p, s) then we would con- 
sider, in place of Eqn.(28), the external field La- 
grangian 

AC = sj^j^sA'' 
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and the corresponding sum rule for the strange 
quark matrix element in Eqn.(37) will have, in its 
left hand side, only the second and fourth terms 
due to modifications of the vacuum and all other 
terms will be set to zero. In the light of our dis- 
cussion of OZI rule, one may then expect that the 
sum rule Eqn.(37) will work better for the isoscalar 
matrix element than for the octet or SU(3) singlet. 

3. The sum rule, Eqn.(37), differs from those of 
ref.[l3[ll also in the coefficient of the third and 
fourth terms due to difference in the calculation of 
Wilson coefficients. This point is elaborated in the 
Appendix. 



IV. EXTERNAL FIELD INDUCED 
CORRELATORS 

We now turn to the computation of the constants F 
and H defined in Eqns. (33-36). First consider the isovec- 
tor case with the current defined in Eqn.(29) and the 
corresponding correlator 

m(0)7,75u(0) - d(0)7,75d(0))|0) (38) 

One can write 

where the coefficient of g^;^, IIi, gets contribution only 
from spin 1"*" states while the coefficient of qptqn^ 112, has 
contributions form both 1+ and 0^ states. 

It is easy to see, using the gauge condition, A^q^ = 
and Eqns. (28), (32) and (34) that 



F(/=l) = -n(=i(q2=0) 



(39) 



This can be easily evaluated using the ward-identity 
which reads 

- Ili=Hq')q'+ni=\q^)q^ 
= {rriu + md){0\uu + dd\0) 

-i{mu + TUdf J d^xe''^-^0\T{dj5uix),uj5dm\0) 

(40) 

Isolating the pion pole in 112(9^) and pseudoscalar corre- 
lation matrix element in the r.h.s. we can write near the 
pion pole 



'Hi=\q')q' 



2„4 



(ml - q^) 
From which we get 



p2 4 

Flml+ , : \. (41) 



[ml - g2) 



(42) 



For the isoscalar case we need to consider the analogue 
Eqn(38) 



n 



/=0 
flu 



- J d'^xe"'-'' 

(0|r(i2(x)7p75u(a;) + d{x)^ ^^zd{x) , 
w(0)7.75"(0) + rf(0)7,75d(0))|0) (43) 



with 



F(/ = 0) = -n(=°(q2 =0) 



(44) 



However, unlike the isovector case the ward identities are 
no longer simple since it involves anomaly Eqn.(lO) and 
Eqn.(40) is replaced, in the right hand side, by consider- 
ably more complex set of terms including the anomaly. 
We shall therefore not use the ward identity to find 

n(=«(o). 

On the other hand n{~°(g^) satisfies a dispersion rela- 
tion. Using the operator product expansion for Ti{^^{q^) 
the value of n{^*'(0) can be found by the QCD sum rule 
method. Analogous procedure of course can be used for 
the isovector case as well. This procedure is well known 
and complete details can be found in ref.|2^. Here we 
shall simply write the final result. Denoting by Lm^ the 
Borel transform [3ll we have from 



ni(g') 

iM2Hi(g2) 



1 [ Qni(s) ds 



subtractions (45) 



7rM2 



s - q- 

5Hi(s) e-"/^^' ds (46) 



where 3Hi denotes imaginary part of Hi. Ean.H46|l is 
usually used to compute the mass of the /i state when 
isoscalar correlator is considered and Ai when isovector 
correlator is considered ,24] . Similarly from 



ni(g2) - Hi(0) _ 1 [ !3Hi(s) ds 
s{s — q^) 



q^ 

we obtain. 



+ subtractions 



L 



M2 



q^ 



ni(o) 

Af2 



(47) 

SHi(s) e-'/^^" ds 
s 

(48) 



which can be used to calculate Hi(0). From the OPE 
expansion for the isovector current correlator Eqn.(38), 
we have 



ni(g2 



1 

47r2 
as 



1 



32 64\ asTT 2 
(49) 



where rh = {m^ + md)/'2.. The above Ean. H49|) has also 
been used in the analysis of r decay (25^. In the right 
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0.018 



0.018 r- 

0.016 - 
0.014 ^ 
0.012 - 

I 0.008 - 
0.006 - 
0.004 - 
0.002 - 
- 

0.8 1 1.2 1.4 1.6 1.8 

FIG. 1: The plot show -n{^^(0) obtained from Eqn.(48) 
in the interval 0.8 < < 1.8(GeV")^. The upper curve is 
obtained for model A, cf. Eqn.(48), and the lower curve is for 
model B. See text. 



hand side of Eaii. (|I5|l and Eqn.l^HJ we use for the Ai, 

the central mass value from experiments and adopt two 
models 

model A : '^Ui{s) ^ 7rh\m\^ 5{s ~ m\j (50) 
model B : ^U,{s) ^ _ ^^^p ^ ^.^^^ (51) 

with r 300 MeV. Following Zyablyuk ^ we use the 
values 

5.'(G^.G;:J =0.5 (GeV)4 (52) 
and for the four quark term 

— a, =3x10-3 (GeV)6 (53) 
y 

Ean. (|53|l takes into account the renormalization correc- 
tions to the 4 quark operator computed by Adam and 
Chetyrkin |2^. We note that the estimate in ean. H53() 
is smaller than the value used in the original work by 
ShifFman, Vainshtein and Zakharov J4J. 
Using the above we find (cf. Fig^ 

ni(0) = -0.0154 (GeV)^ (model A) (54) 

ni(0) = -0.0138 {GeV)^ (model B) (55) 

which is to be compared with the ward identity value 
Eqn.(40) with = 0.017 (GeV)^. We like to stress 
that the above calculation does not use ward identity 
and PCAC. 

Consider now the isoscalar case 

w7.75"(0) + d7.75d(0)}|0) (56) 

with 



0.016 ~ 

0.014 - 
0.012 - 

c 

' 0.008 - 
0.006 - 
0.004 - 
0.002 - 

I ^ ^ ^ ^ 

0.8 1 1.2 1.4 1.6 1.8 

FIG. 2: n{=°(0) obtained from Eqn.(48) using /i(1285) and 
ignoring its width as in model A. 



Before proceeding with the calculations we note 
the interesting similarities between the vector states 
/9(770), cj(780), 0(1020) and the axial vector states 
Ai(1235),/i(1285),/i(1420). Both p and Ax are broad 
resonances and are nearly degenerate with their 7 = 
partners lo and /i respectively. Morever the decays of 
(/)(1020), /i(1420) are dominated by strange mesons. To 
evaluate n(=°(0) we can proceed in an analogous manner 
as for the isovector case above. First we note that in the 
physical side or the r.h.s. of the sum rules, Eqns.(46) and 
(48), /i(1285) will replace Ai. The higher mass state 
/i(1420) is effectively included in the sum over higher 
mass states which in the usual QCD sumrule approach 
are represented by the quark loop contributions with an 
effective threshold W'^ . We should stress here that the 
details of the decay modes of /i(1285) and /i(1420) are 
irrelevant and do not enter the sum rules. 

Now turning to the OPE it is clear, from Eqn.(38) and 
Eqn.(43), that as long as quark anihilation diagrams are 
neglected the OPE for 1=1 and / = currents will 
be identical. Consider the first term in Eqn.(49). This 
arises from the single quark loop and is the same for the 
isovector and isoscalar case. However in the isoscalar 
case we must include corrections that can arise from a 
three loop diagram in which the initial quark current 
loop annihilates to a two gluon intermediate followed by 
materialisation to the final current quark loop. These 
diagrams are expected to contribute with coefficients like 
a^/TT^ and for that reason we neglect them here. 

The calculation of n(^'^(0) is then straight forward and 
we obtain (cf. Fig [21) 

n(=°(0) = -0.0152 {GeVf (57) 

In the analysis of the sumrule Eqn.(37) we shall therefore 
use F = 0.0152 {GeVY as the central value and study 
the effect of variation around it. 

We now turn to the determination of the dimension 
5 correlator or the constant H. For the isovector case 
this has been computed by Novikov et.al. i27(|, using the 
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sumrule for 112 in 

= -Ili{q^)g^, + n2{q^)q^q, (58) 
defining 

{0\gJ^G^^^jpu{x)\7T) = -S^F^k^. (59) 

They obtained — 0.21{GeV)'^. However they had used 
a somewhat high value for the four quark correlator. This 
has been reanalysed by lofFe and Oganesian using 
more upto date values of various paramters. They obtain 
(5^ = 0.16{GeV)^ . We have independently reanalysed the 
sum rule of Novikov et.al. 27] and agree with loffe and 
Oganesian 15]. From Eqn(58) and Eqn(59) one gets the 
value H{I = 1) = 6^F^. 

To find the value H in the isoscalar case we proceed as 
follows. loffe and Khodjamirian fisl] have computed 

(0|5sX!^Y'^/3"'^'3'?(^)|0)Is;7(3) szriglet = (60) 
9 

and find ho « 3x 10~^ (GeV)^. Since the matrix elements 
in Eqns.(59) and (60) are quite small we can use SU(3) 
flavour symmetry and write 

2 A" ~ 

= 3(0|5s9yG'"^7<^'7|0)|s;7(3) singlet 

+ l{0\g,q^G:^^,q\0)\o,tet (61) 

From which we obtain 

H = ho + y^F^ 

« 1.14 X 10"^ (GeV)^ (62) 

We shall use this as the central value in the our analy- 
sis of the sumrule Eqn.(37). The above term is far less 
significant in the determinations of G in Ean. (|37|l than 
the dimension 3 correlator constant F calculated in the 
earlier part of this section. 

V. DETERMINATION OF Gu + Gd FROM THE 
SUMRULE 

We can use the values of F and H determined in the 
previous section in Eqn.(37) to find Gjj + Gd- But first 
we need to determine which can be obtained from 
loffe's sumrule for the nuclcon mass. 

E2 bM^ Eo ia^L'^^^ a^ml 
L4/9 + 4L4/9 3 3M2 

= A^ e-™«/*^' (63) 



1 1 ^ ^ ^ 1 

0.8 1 1.2 1.4 1.6 1.8 

FIG. 3: A^ obtained from Eqn.(63) is shown in the Borel 
mass interval 0.8 < A/^ < 1.8 {GeVf. 



0.2 




0.17 I ^ ^ ^ ^ ^ ^ ^ ^ ^ 1 

0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 

m2 



FIG. 4: G + AM'^ obtained from Eqn.(37) is shown in the 
Borel mass interval 0.8 < < 1.8 (GeV)^ A straight line 
fit, dotted line gives G = 0.22 and A = -0.026. 



Experimentally we have Aqcd = 247 MeV [30], which 
corresponds to as{l GeV) — 0.5. To find Xj^ we shall 
use the experimental value m^r = 0.938 GeV and we fix 
= 2.22 by looking for the best fit in the least square 
sense in the Borel mass interval 0.8 GeV^ < < 1.8 
GeV2, we find (cf. FigEJ 

A^ 1.975 GeV^ 

This can now be used in the sumrule Eqn.(37), along 
with the values F = 0.0152 (GeV)'^ and H = 1.14 x 
lO^'^ (GeV)^. We find fitting the sum rule with W'^ — 
2.22 and the same interval for as used in Eqn.(63) 
(cf. Figg 

G = Gu + Gd = 0.22 and A = -0.026 {GeVy^ (64) 

Since our calculation of F and H in the previous section 
are not exact we have varied them by 20%. An numerical 
increase of F by 20% yields a value G = 0.34 and A = 
— 0.011(6*6^)"^, while a decrease by 20% gives a value 
G = 0.10 and A = -0.041(Gey)-^ Compared to this a 
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change of value of H by 20% barely changes the result 
by 2%. 



VI. CONCLUSION 

To evaluate the linear combination ^Gu + |G_d + ^Gs 
occuring in the Bjorken sum rule we can proceed as fol- 
lows. We write 

Gsj ~ -^{Gu ~ Gd) + -^{Gu + Gd) 



1 

18 



(G(7 + Gd — "^Gs 



(65) 



First term is known from neutron /3 decay and we have 

G[/-Gd = 1.267 (66) 

The last term in the Eqn.(59) is also known from hypcron 
/? decay assuming the validity of flavour SU(3) symmetry. 
We have [1 



Gu + Gd- 2Gs ^ 0.585 



(67) 



We can now use the value of Gu + Gd determined in 
SecE|Eqn.(64) to get 

G'Bj(Ai^ = 1 GeV'^) = 0.32 (68) 

to be compared with experimental value 

Gbj (^2 ^ 5 Gey2) ~ 0.28, (69) 

which increases slightly when fj,^ is decreased from 5 
GeV^ to 1 GeV^. Taking into account the dependence 
of the singlet matrix element [28L l29l | the experimental 
number, Eqn.(69), increases by a few percent to 0.29 at 
fj,^ = lGel/2. 

In arriving at Eqn.(68) we have used SU(3) flavour 
symmetry via Eqn.(67) in Eqn.(65). Returning to 
Cabibbo theory, all the octet current matrix elements be- 
tween the octet of baryon states are expressible in terms 
of two irreducible matrix elements F and D, and one has 



Gr 



Gu -Gd = F + D 
Gd ~2Gs ^3F - D 



(70) 
(71) 



In QCD sumrule calculations of the various octet current 
matrix elements between baryon states one first considers 
the limit of massless quarks which is of course automat- 
ically SU(3) symmetric. The constants F and D can be 
determined in a variety of ways. For example the ma- 
trix elements of the isovector current u^fij^u — d'y^j^d 
between nucleon states gives F + D while, S ^ S is pro- 
portional to F, E — > A is proprtional to D and S — > S to 
D — F. Details can be found in ref.[l,0|. In ref.Q, using 
a different Lorentz structure invariant for the propaga- 
tor in Eq.(26), the relation IF w hD was found which 
is not far from experiment. To incorporate the effect 
of quark masses one expands in quark mass and retains 



terms linear in them and also take into account the dif- 
ference between the strange quark condensate and up or 
down quark condensate. We first note that since up and 
down quark masses are negligible, the isovector current 
matrix element in the nucleon is unaffected. In a recent 
update of the earlier calculation in ref. jSj , loffe and Oga- 
niesen find 

Gu-Gd = 1.37 ±0.01 

Finding D, F and D — F from the hyperon matrix el- 
ements, when quark masses are included, is sensitively 
dependent on the Borel mass region as was found by 
authors of 9J. For this reason it is difficult to decide 
how accurately Cabbibo theory is satisfied by using QCD 
sum rules. We have therefore relied on experiment which 
is reflected in the use of Eqn.(67) to compute Gsj in 
Eqn.(65). Nevertheless, it is worth emphasising that 
QCD sumrules provide, though not precise, a quantita- 
tive explanation of Gbj, without any arbitrary parame- 
ter and use only the value of vacuum condensates already 
obtained through other hadron properties. 



APPENDIX A: WILSON COEFFICIENTS 

There are two differences between the coefficients of 
the term in Eqn.(37), namely, the coefficient of the 
gluon condensate term, 6, and the coefficient of the ex- 
ternal field induced correlator, H, between ref. [ol Ull 
on the one hand and ref. jlOl [l5l| on the other. Full detail 
of the calculation can be found in the ref.P|. These cal- 
cualtions have since been verified again by the original 
authors of ref.j^ themselves and authors of ref.^^ and 
ref. ^2 ■ In ref. Q the nucleon current correlator, Eqn(26) , 
is calculated for a generic A£ 

A£ = (5'„U7m75M + 9ddl^Ll5d)A'^ 

with gu and gd as arbitray constants. According to Table 
1 in ref. 0, the contribution of Fig (5) and (6) in ref.0 
have coefficient —{2gd + 10g„/3) and 2(g„ -I- gd)- In the 
computation of Ga, Qu = —gd so that Fig (6) gives zero 
and there is complete agreement between ref. ^ and 
However for the isoscalar current, for which g„ = 
gd = I, the results are different. Ref. [ajll3 have 



~(2.gd+yg«)+2(.g„-Hgrf) = --.g„ 
Reversing the sign of Fig (6) gives 



(.gu-ffd) (Al) 



10 , N 28 

-(2gd + yffu) - 2{gu + gd) = -yfu [gu = gd) 

(A2) 

which is r esult of ref. [13. Thus ref. 9] and HO] agree 
numerically in the case gu — —gd but different in the 
isoscalar case. In Eqn.(37) the coeffcient used corre- 
sponds to Ean. l|Al() . 

We now turn to the coefficient of the gluon condensate 
{glG'tyG'ti) ■ Again details of the calculations are given 
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in rcf. 9] in their Table 1 which uses their cqn(2.12) and 
Fig. 2 which leads to coefficient —g^/A. On the other 
hand ref. llOt s eem to have gd/i for this coefhcient. Thus 
ref.0 and [lOf agree numerically in the case g„ = —gd but 
have opposite signs in the isoscalar case. Interestingly 
enough if one tries to obtain the Wilson coefficients for 
T{'r]{x)f]{0)}\A^ by using a chiral rotation of the quark 
fields given by 

u — > e'f" -^-^ 
d — , e'S" -^-^ '^"d 

in the propagator (0|r{?7(a;)77(0)}|0) in absence of the 
external field, which occurs in the mass sum rule and 
assuming the validity of chiral invariance of the vacuum 



to obtain {0\T{ri{x)fj{0)}\AjO) one obtains gd/4: for the 
coefficient of (gsGI^^GIj^^,). On the other hand it is well 
recognised that presence of non-perturbative gluon fileds 
in the vacuum state leads to breakdown of chiral symme- 
try, which would seem to invalidate such a calculation. 
Again in Eqn.(37) we have used the result — (7„/4 found 
in ref.Q and ref.|0|. Fortuitiously enough these two 
differences in the coefficients of third and fourth terms 
in Eqn.(37) between ref.0 and ^3 numerically tend to 
compensate. Also we have seen the dimension three term 
characterised by F which occures with M"^ in Eqn.(37), 
is far more significant than the Af^ term in obtaining 
Gjj + Gd- Our main point has been that F can be com- 
puted from the spectrum of axial mesons and we get a 
sensible answer for Gjj + Gd and therefore G bj ■ 
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